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Abstract 

In this paper we classify and derive closed formulas for geometric elements of quadrics in rational Bezier triangular form 
(such as the center, the conic at infinity, the vertex and the axis of paraboloids and the principal planes), using just the 
control vertices and the weights for the quadric patch. The results are extended also to quadric tensor product patches. 
Our results rely on using techniques from projective algebraic geometry to find suitable bilinear forms for the quadric 
in a coordinate-free fashion, considering a pencil of quadrics that are tangent to the given quadric along a conic. Most 
of the information about the quadric is encoded in one coefficient, involving the weights of the patch, which allows us to 
tell apart oval from ruled quadrics. This coefficient is also relevant to determine the affine type of the quadric. Spheres 
and quadrics of revolution are characterised within this framework. 

Keywords: Algebraic projective geometry, rational Bezier patches, quadrics, Steiner surfaces. 


1. Introduction 


One of the reasons for extending the framework of sur¬ 
faces used in CAD from polynomial or piecewise polyno¬ 
mial surfaces to rational ones is the inclusion of quadrics 
in an exact fashion, since quadric surfaces such as spheres, 
cones, cylinders, paraboloids are commonly used in engi¬ 
neering and in architecture [ij. 

There are mainly two ways of implementing surfaces in 
CAD: tensor product patches and Bezier triangles 0. The 
former ones are the most common, but the latter have 
applications in animation and finite element theory due to 
the flexibility of triangles, compared to quadrilaterals, for 
constructing surfaces and adapting to different topologies 
without producing singularities. 

Whereas rational quadratic curves are conics, we have a 
different situation when we go to surfaces. On one hand, 
not every quadric triangular patch can be represented as 
a rational quadratic Bezier triangle. On the other hand, 
rational quadratic Bezier triangles are in general quartic 
surfaces H and only in some cases quadric patches are 
obtained y. 

Quadric triangular patches have been studied from 
many points of view. [3 describes a method for construct¬ 
ing quadric patches grounded on a quadratic Bezier con¬ 
trol polyhedron. @ and 0 use an algebraic Mproach to 
construct curves and surfaces on quadrics. (8| classifies 
quadratically parametrised surfaces and studies their ge¬ 
ometry. 0 provides a thorough classification of rational 
quadratic triangles. provides an algorithm for check¬ 
ing whether a rational Bezier triangle belongs to a quadric 
and classifying it. This solution is shown in ll| to be one 


of fifteen possibilities, with different numerical condition¬ 
ing. In [12| the algorithm is used to obtain the axes of 
the quadric, [l^ provides a tool for constructing ratio¬ 
nal quadratic patches on non-degenerate quadrics. Three 
corner points and three weights are used as shape param¬ 
eters in 1^ to design quadric surface patches. In 15| the 
shape parameters are three points and the normal vectors 
at them. 


In this paper we draw geometric information from ratio¬ 
nal Bezier quadric triangles and tensor product patches. 
We classify them and calculate their geometric elements 
in closed form, using just the control net and the weights 
of the patch, as it is done for conics in [l^. One of the 
reasons for doing this is that closed formulas for geometric 
elements of a quadric can be of great help for designing 
with them. 


With this goal in mind, we derive bilinear forms for 
the quadrics, both in point and tangential form in a 
coordinate-free fashion, using techniques borrowed from 
algebraic projective geometry, such as the use of pencils of 
quadrics, which have also been used in [l3)[H)[l3|- Using 
linear forms with clear geometric meaning (tangent planes, 
planes containing boundary conics...) instead of, for in¬ 
stance, cartesian coordinates, allows us to find a closed 
form for the implicit equation of these quadrics in terms of 
their weights, their tangent planes and the plane spanned 
by the corner vertices of the patch. The major originality 
and advantage of our formulation is the use of only geo¬ 
metric information about the patch, avoiding the use of 
coordinate expressions as an intermediate step. In fact, 
most of the relevant geometric information of the surface 
which is necessary for its classification is comprised in a 
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single coefficient, the parameter of the pencil of quadrics, 
which depends on the weights of the patch. 

The coordinate-free bilinear forms for the quadrics en¬ 
able the derivation of closed formulas for several geomet¬ 
ric elements (center of the quadrics, bilinear forms for the 
conic at infinity, diametral planes) calculated in terms of 
weights and vertices. To our knowledge, such closed for¬ 
mulas have not been produced before. Similarly, linear 
forms for principal planes are obtained up to solving a 
cubic equation. The degeneracy of the solutions of this 
equation allows us to identify quadrics of revolution and 
spheres. 

Finally, the previous results for triangular quadric 
patches can be extended to rational biquadratic quadric 
patches using the tangent planes at three corners of the 
boundary of the tensor product patch. 

This paper is organised as follows: We revisit ratio¬ 
nal Bezier triangular patches and the characterisation of 
quadric patches in Section 2 . In Section 3 we construct the 
pencil of quadrics, in point and in tangential form, which 
are tangent to our quadric at the conic through the three 
corner vertices of the patch and fix the only free parameter 
in terms of the weights of the quadric. This provides us bi¬ 
linear forms for the quadrics and hence implicit equations. 
A coordinate-free expression for the center of the quadric 
is obtained in Section 4 as a barycentric combination of the 
corner vertices of the patch and the intersection of their 
tangent planes. It is used to tell paraboloids from centered 
quadrics. In Section 5 we calculate a bilinear form for the 
conic at infinity, which is useful for determining the affine 
type of the quadric. Section 6 is devoted to degenerate 
quadrics. With this information we provide in Section 7 a 
way to classify quadric patches using the signature of the 
bilinear form, the center and the bounding conic arcs. We 
introduce the scalar product in Section 8 in order to calcu¬ 
late Euclidean geometric elements such as principal planes 
and axes of quadrics and vertices of paraboloids. This al¬ 
lows us to characterise spheres and quadrics of revolution. 
In Section 9 we show several examples of application of 
our results. We show how to extend our results to quadric 
tensor product patches in Section 10 . A final section of 
conclusions is included. 


such that u + V w = 1 , u,v,w € [ 0 , 1 ], where the coeffi¬ 
cients Cijk are the control points for the surface and uiijk 
are their corrresponding weights. 

The surface patch is bounded by three rational Bezier 
curves of degree n, which are obtained by fixing u = 0, 
V = 0 , w = 0 . Their respective control points and weights 
are respectively the ones with index i = 0, j=0, fc = 0. 

We are interested in the special case of quadratic sur¬ 
faces. In this case the control net and the matrix of weights 
are 
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and the surface patch is completely determined by the 
three boundary conic curves: 

The conic at m = 0 has control points {coo2 ,cqh, co2o} 
and weights {uJoo2,<^oii,<^02o}, the one at z; = 0 has con¬ 
trol points {coo2,cioi,C2oo} and weights {a;oo2, wioi,a;2oo}, 
and the one at w = 0 has control points {co20) cuo, C200} 
and weights {wo20: ^^iiO: ^^200}- We just consider the case 
of non-degenerate boundary curves. We name the planes 
where such conics are located respectively as u,v,w. 

Quadratic rational Bezier triangular surfaces comprise 
quadrics as a subcase, but in general they are quartic 
surfaces named Steiner surfaces Q. On the contrary, 
not every quadric patch bounded by three arcs can be 
parametrised as a quadratic rational Bezier triangle. Ac¬ 
cording to 0: 

• If the Steiner surface is a non-degenerate quadric, the 
three conic boundary curves meet at a point S and 
their respective tangent vectors at S define a plane. 

• If the three conic boundary curves meet at a point 
S and their respective tangent vectors at S define a 
plane, the Steiner surface is a quadric patch. 

• If there are three alligned points, one on each bound¬ 
ary conic curve, and the tangent plane to the surface 
is the same for all three points, the Steiner surface is 
a degenerate quadric patch. 


2. Quadratic rational Bezier triangular patches 


Rational Bezier triangular patches of degree n are de¬ 
fined using trivariate Bernstein polynomials of degree n in 
parameters u, v, w. 




nl 
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i + j + k = n, 


. , i+j+k=n , . 

c{u,v,w) = -—-—-—, ( 1 ) 


The last case just implies that the quadric (cylinder or 
cone) is ruled and the three points lie on the same ruling. 

We are considering Steiner surfaces with such a point 
S. This case comprises non-degenerate quadrics. Since 
S belongs to the three boundary conics, there are values 
Us, Vs,Ws of the parameters such that 
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where we have introduced quadratic Bernstein polynomi¬ 
als, 

However, since the set of weights of rational Bezier 
curves is unique up to Mdbius transformations [l^ of the 
parameter, we can use this freedom to set Ms, Vg, Wg at in¬ 
finity. It is a simple exercise to check that this is possible 
if and only if the tangent vectors at S are coplanar [I^ . 
Hence, this choice is not a restriction. 

From now on we assume that the set of weights for the 
surface fulfills this condition and therefore we may write 
S in three different forms as 

^ _ W002C002 — 2 a;oiiCoii -I- W020C020 

UJu 

_ W002C002 — 2 a;ioiCioi + W200C200 

UJv 

_ W200C200 ~ 2 a;iioCiio + W020C020 

— I I'JI 

UJyj 

where uiu, lOvi^w guarantee that these are barycentric com¬ 
binations, 

Wu = 02002 — 2 o;oil + 02020: 02„ = 02002 — 2 o 2 ioi + 02200: 
Ulw = 02200 ~ 2o2iio + 02020: 

if S' is a proper point. If S is point at infinity, the coeffi¬ 
cients 0Ju,uJv, ujuj are chosen in order to have the same vec¬ 
tor as representantive for S on the three boundary curves. 


3. Pencils of Steiner quadrics 


In order to obtain a bilinear form for the quadric, we 
describe the pencil of quadrics which are tangent to our 
quadric at the conic on the plane t defined by the points 
at the corners of the surface patch, P = coo 2 : Q = C 020 : 
R = C 2 (ro. This pencil of quadrics has been used also in 

BEQ- 

To this aim, we just need two quadrics belonging to the 
pencil [13, which can be degenerate. 


Figure 2: Conic circumscribed by a triangle 


Clio 



Figure 3: Tangent planes p, q, r 


For instance, we can take the cone C tangent to our 
quadric along the conic on t and the double plane t (see 
Fig.p. Hence, the bilinear form for this pencil of quadrics 
is just C + using the same letter for a surface and its 
form. Finally, we can determine our quadric just using the 
existence of point S, 

C{S) + \t{Sf = 0. 

We just have then to provide a bilinear form for the cone 
C. Since A = —C{S)/t{S)^, we notice that the sign of A 
depends on whether S lies in or out of the cone C. 

It will be useful to start with the bilinear form for the 
conic on t (see Fig. [ 21 ). In this figure u,v,w are the straight 
lines which are respectively the intersections of the planes 
u,v,w (where the boundary conics lie) with the plane t. 

We use the same letter for a point (uppercase) and its 
polar line (lower case). If a point A lies on the conic, 
its polar line is the tangent a to the conic at such point. 
Hence, p, q, r are respectively the tangent lines to the conic 
at P, Q, R (see Fig.|31). We shall use the same letters p, q, r 
also for the tangent planes at P,Q, R. 

On the other hand, the polar line to a point A not on 
the conic is the line a linking the tangency points of the 
tangent lines to the conic drawn from A. Hence, u, v, 
w are the polar lines to the vertices of the circumscribed 
triangle U, V, W. 

The pencil of conics circumscribed by the triangle de¬ 
fined by the lines p, q, r is easily described in tangential 
form. The bilinear form for this pencil is just 

aUV + PUW + 7 HIF, 
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with coefficients a, /3, 7 . It is easy to check that this form 
vanishes on p, q, r, since 

p{U) = 0 = p{V), q{U) = 0 = q{W), r{V) = 0 = r{W). 

But, for our purposes, we require the point bilinear form 
for the conic. Referred to the lines r, q, p, which form a 
dual reference to the one formed by U, V, W, provided 
that their linear forms satisfy 

piW) = l, qiV) = l, r(C/) = l, (4) 

the matrix of this point bilinear form is the inverse of the 
one for the tangential bilinear form, 

7 ^ —/37 —a'y \ 

-h -a/3 
—07 —a/3 a^ J 

and hence such a bilinear form is 



2af3j 


a^p^ — 2al3pq — 2a'ypr + — 2/3 jqr + 


' 2„2 


„2„2 


In order to determine the coefficients a, /3, 7 , it is useful 
to take into account that provides relations between 
values of the linear forms for p, q, r, 

= P{S) = 


W, 
W 2 OO 
UJyj 
^002 


q{R) = q{S) = 
r{P) = r{S) = 


LUy UJu 

which can be used to write U, V, W in terms of P, Q, R, 
W 002 Wu)P + 0J020^vQ — I^200^uR 


u = 


V = 


fiu 

W 002 Wii)P — i^020^vQ 


^200‘^uR 


fly 

— W 002 Wu)P + ^020<^vQ - 


- ^20Q^uR 


n 


w 


W = 

or, conversely, 

^ iljjU + ^lyV 

2 woo 2 Wu, ’ 2 a;o 20 <^ 

where the normalisation denominators are 


( 5 ) 


Q.wW + fljjU Q,vV + Q-wW 
; 0 : 7^ ; -ft — , 

ZUJ200i^u 


flu = <Uoo2<-^w + 0J020^v — W200Wuj 
fly = UJ002<Z’w — (^020(^v + W200^^ti: 
flw = — W 002 WU, + a;o2oWt, + uj2oo^u, (6) 

provided that U, V, W are proper points. If one of these 
points goes to infinity, the respective denominator is taken 
to be one. The interpretation of these coefficients as 
weights for U, V, W on the conic on t is discussed in Sec¬ 
tion 6 . 

Finally, the normalisation condition ([4]) allows us to fix 
the linear forms for p, q, r since 

P(W) = g(V) = 

,.([/) = ?^r{S). (7) 

llu 


Since we require that points P, Q, R lie on the conic, 
l3q{P) = 7 r(P), ap{Q) = 7 r(Q), ap{R) = l3q{R), 


we identify the coefficients of the bilinear form. 


n 


w 


/3 = 


fl\ 


7 = 


fir 


up to a common factor. 

Hence, a bilinear form for the conic on t is 


p^ q^ 2pq 2pr 2qr 

fl^ fly fffj fl\yfly flypflij flyfljj 

and, in tangential form. 


fluflyUV -l- fluflwUW flyflw^W. (9) 


The terms of the bilinear form can be factored as 



and thereby it has signature (- 1 - H-) regardless of the 

values of the coefficients. 

We move now back from plane to space: If p, g, r are the 
tangent planes to the quadric at P, Q, P, the bilinear form 
m describes a degenerate quadric of signature (- 1 - H— 0 ) 
which contains the conic on t and it is tangent to our 
quadric along it. Since the rank of the bilinear form is 
three, it is a cone (or a cylinder, if the intersection point 
of p, q, r goes to infinity), since these are the only non¬ 
plane degenerate quadrics. Hence, it is the cone we are 
looking for. 

The bilinear form for the pencil of quadrics referred to 
the planes r, q,p, t, 


_L ^ I 2 

fl^ fly fl^ flyyfly flyyflu flyflu ' 

produces a bilinear form for our quadric if A = 

-C{s)/tisr. 

We see that the sign of A determines the signature of the 
bilinear form, defined as the difference between the number 
of its positive and negative eigenvalues. This is relevant 

since ruled quadrics have null signature, (-|- H-) and 

oval quadrics have signature two (- 1 - H-h). This is useful 

to classify the quadric. 

In tangential form referred to the points U, V, IT, T, the 
bilinear form for the quadric is 

rp2 

Q.jjQ.yUY + Q:JjQ:\yUW + Q.yQ:wYW -—, 

A 

if t{T) = 1, where T is the intersection point of the planes 
p, q, r, that is, the vertex of the cone C. 

With this information we can write T as 


T = 


+ 


^ 002 ^u;(^u T ^w'/R 3~ U^020^7j(^7i “f ^w')Q 


euflu^u eyflyiOy I eypflyp(jJpp 2eQLyipiOyLO' 

^200^'u( 3~ UJy UJyy'jR 2UJuUJyUJyjS 


( 10 ) 
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4. Center of a quadric 


or alternatively, 


5 = 


n 
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■U ■ 


0. 
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-F + 
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W- 


2 esw y^tUyUJ-Kj ClJ^UUJy C\/^Y^V C\Y^\YUJy 

‘2UJq>UJyUJ-tn 


T, 


where takes the value one if the corresponding point A 
is proper or zero if it is a point at infinity. 

Now we can compute A and determine the quadric: 


The bilinear form for the quadric provides a way to ob¬ 
tain its center Z as the pole of the plane z at infinity. Since 
the plane at infinity is formed by vectors, we may write its 
elements as barycentric combinations 

all + bV + cW + dT, a + b + c + d = 0. 


Theorem 1. A rational triangular quadratic patch with 
control net {coo 2 ,coii,co 2 o,cioi,cuo,C 200 } and weights 
{woo 2 ,Won,^ 0207 ^ 101 ,^ 1107 ^ 200 }, such that the three 
boundary conics meet at a point S, which is written as 
in m, is a quadric with bilinear form 


We have to consider the possibility of having any of the 
points of the reference at infinity. In such case, the null 
sum is restricted to proper points, 

all + bV + cW + dT, eua + eyb + ewc + CTd = 0.(14) 


Q 


p^ q^ 2pq 2pr 2qr 2 

fly fffj flwflv flw^U flv^U ’ 


where t is the linear form of the plane containing 
C002, C0207 C2007 P is the linear form of the plane containing 
C 0027 Con, cioi, q is the linear form of the plane contain¬ 
ing C020, Con, Clio and r is the linear form of the plane 
containing C200, cioi, cno which satisfy 


p(W) = l, qiV) = l, riU) = l, t{T) = l, 


at the points U (intersection of the planes p,q,t), V (in¬ 
tersection of the planes p,r,t) and W (intersection of the 
planes q,r,t) given by (0) and the point T (intersection of 
the planes p, q,r), given by m- 

The coefficients are given in except for X, 

^ 4 ” 4 ” 2 UJyUJyy 2 UJyjUJy 

(2esuJuajyUjyj — eufltjuiu — — ewflwOJw)'^ 

and the tangential bilinear form for the quadric is 

rp2 

Q = Q.i/Q.yUV + Q.ijfl\YUW + , (13) 

A 

Furthermore, for A > 0 the quadric is oval and for A < 0 
the quadric is ruled. 


The expression we have obtained for the bilinear form 
of the quadric has the advantage of encoding most of the 
information about the surface in just the coefficient A for 

e. 

This result provides a procedure for computing a bi¬ 
linear form, and hence the implicit equation, for a non¬ 
degenerate Steiner quadric patch in a coordinate-free fash¬ 
ion using just the vertices of the control net and their re¬ 
spective weights: 

1 . Obtain S as intersection of the planes u,v,w. 

2. Compute the normalised linear forms for the planes 
t,p,q,r. 

3. Obtain an equivalent list of weights fulfilling ([3]). 

4. Use Theorem [T] to obtain the bilinear form Q for the 
quadric patch. 

5. The implicit equation for the quadric patch is then 
Q{X,X)=0. 


Hence, a linear form for the plane at infinity in this 
reference is just 


2 = ewP + evq + ec/c + ert. 


(15) 


The pole of this plane, 


Q{z) = ewflw{fd,uU -\- UyU) -|- -\- fluU) 

-\- eijflu{flvV -\-flwW) -^2^7 (16) 

A 

can be written in a simpler way in terms of P, Q, R, T 
in order to produce an expression for the center of the 
quadric. 




ewflwa!waJoQ2P + evflv^y<^02oQ + eu^ua>uaJ20oR 


Ct 

Xflz 


T, 


flz 


(17) 


where the denominator flz is one if the center is a point 
at infinity or, if it is a proper point. 


flz 


flwa>wa!oo2 + flva!vUJo2o 4- flua!uUJ2oo 


er 
A ’ 


(18) 


and we have introduced for simplicity, 


11(7 = cc/Hc/, Hy = Eylly, Hw = ewHw- (19) 

Since the center of a paraboloid is a point at infinity, we 
have a simple characterisation: 

Corollary 1. A rational triangular quadratic patch with 
control net {0002,coii,C0207cioi,ciio,C200} and weights 
{u;oo 2 ,woii,wo 207 Wioi,wiio,W 2 oo}, such that the three 
boundary conics meet at a point S, which is written as 
in m, is a paraboloid if the quadric is non-degenerate and 

~ ~ Cx 

llvyWi(,a;oo2 4- Hya;„a;o20 4- fluUJuUJ200 —~ 

where the coefficients are given in W, IW and m- 
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5. The conic at infinity 

The conic at infinity is the intersection of the quadric 
with the plane 2 at infinity and it is formed by its asymp¬ 
totic directions. It is useful for classifying quadric patches, 
as it is done in [T^. Since points X on z satisfy 

z{X) = ewp{X) + evq{X) + eur{X) + eTt{X) = 0 , 
we can use as bilinear form for the conic at infinity on z 
^ p2 ^2 ^2 2 pq 2 pr 2 qr 

fly fly flv^fly flu'flf/ flyflf/ 

+ A(eu/p + cy? + (20) 

except when T is a point at infinity. 

In order to draw information about the conic at infinity, 
we may factor its bilinear form, 


an ellipsoid then, since it does not intersect the plane at 
infinity. 

On the other hand, if the determinant m is negative, 
the signature of the bilinear form is (-1- H—) and we have 
a proper conic. The quadric is a hyperboloid in this case: 

Corollary 2. A rational triangular quadratic patch with 
control net {coo2,coii,co2o,cioi,cuo,C200} and weights 
{woo2,Won,W020,wioi,W ho,W200}; such that the three 
boundary conics meet at a point S, which is written as in 
m and with positive X, given by Qgl) . is an oval quadric 
and: 

• An ellipsoid ifdetZ > 0 . 

• An elliptic paraboloid ifdetZ = 0 . 

• A two-sheeted hyperboloid ifdetZ < 0 . 


Z 


(O^A -I- l)(Ti -I- -2 , . , ^2 
fl^A-b 1 

-|- tlyClw + Oy/Oj/ — 



where the expression for det Z is given in m- 

6. Conics on planes t, u, v, w and degenerate 
quadrics 


where we have introduced three linear forms, cri, <72, era, in 
order to diagonalise Z, 


CTi := 


p flyflvvA — 1 q flijflwX — 1 r 
Xlw O^A -1-1 fly fiwA -1-1 Xlu 




(72 := (fly + fliy) 


n 


y 


-b 


f2[/fly + flyfly/ -b fly/flfj — fl^ — 2/A r 


fly -b fly/ 


Or/’ 


0-3 := 


Oy -b fly/ 0(7 


The case of T at infinity is simpler, as points at infinity 
satisfy 


z{X) = ewp{X) + evq{X) + eur{X) = 0, 


but it can be handled similarly. 

Combining both cases, we obtain a general expression 
for the determinant of Z in this reference, 

det Z = A ^OryOy -b OyOy/ -b Oy/0(7^ — Ct- ( 21 ) 

The conic at infinity of a paraboloid is degenerate. 
Hence, detZ vanishes for these quadrics. This condition 
is equivalent to the one obtained in the previous section. 

The classification of the conic at infinity allows us to 
finish the classification of quadrics. Since paraboloids are 
non-centered, one-sheeted hyperboloids are ruled and cen¬ 
tered, we just have to tell ellipsoids from two-sheeted hy¬ 
perboloids, since they are both oval and centered: 

Since A is positive in this case, we notice that if the 
determinant (l?T]) is positive, the signature of the bilinear 
form is (-b-b-b) and hence the conic is imaginary. We have 


In order to draw more information about the surface 
patch, we take a look at the conic arcs on planes t,u,y, w. 

A conic with weights a/o,wi,a/2 can be classified [l8| us¬ 
ing the canonical weight uj = wi!y/iZfuJf,'. 

If w < 1 , it is an ellipse; if w = 1 , it is a parabola and if 
w > 1, it is a hyperbola. 

The canonical weights for the conics on planes u,v,w 
are respectively 

Won wioi Who 

\/W002W020 ’ \/Wo02W200 ’ \/W020W200 

A set of weights for an arc of the conic on t is readily ob¬ 
tained. For instance, we can use P, C/, Q as control polygon 
and use the same kind of construction of (| 31 ) to assign an 
infinite parameter to the point R, so that its coordinates 
referred to P,U,Q, 

W002Wu,P — H/zt/ -b UJq20U!vQ 

R = -, 

W200W11 


provides us a set of weights woo2Wu,, H(7/2, 0/020Wu for the 
conic arc. 

Similarly, we obtain weights 0/002Wm, Hy/ 2 , 0/200W11 for 
the arc with control polygon P,V,R and o/o2ow„, Hy// 2 , 
W2ooWu for the arc with control polygon Q, W, R. Hence 
the canonical weights for these arcs are 


Hr 




WPC/Q = 


2 v^W002W020W)A^ 
a^QWR = 


l^PVR = 


H 


w 


2-v/W020W200WiiW„ 


2-\/W002W200WuWu, ’ 

( 22 ) 


We can use either of these to classify the conic on t. 
Furthermore, this result furnishes an interpretation of 
flu, fly, flw as weights for the points U, V, W if 0/002w^,, 
wo2oWu and o/2ooWu are respectively the weights for P, Q, R. 
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We have calculated this set of weights resorting to the 
point S', but it is clear that it can be obtained indepen¬ 
dently from the barycentric combinations 

p _ ^IqQ — riwW + vlrR _ ripP — -i- ripi? 

Oq — rirV “I" ^-R ’ “f ^-R 

n _ OrP - flpU + OqQ , 

rip-nu + ^iQ ’ ^ ' 

up to a multiplicative factor. 

This is useful for degenerate quadrics (cones and cylin¬ 
ders), since for their triangular quadric patches the bound¬ 
ary conics do not meet in general at a point S, but we can 
still use ([U as the bilinear form for the tangent cone to 
the quadric along the conic on t: 

Theorem 2. A rational triangular quadratic patch with 
control net {coo2,coii,C020,cioi,cno,C200} and weights 
{woo2,woii, ^0207^101,^1107^200}, on a degenerate quadric 
has a bilinear form with coefficients ffc/, ily, flw satisfying 

m, 

n _ I I (241 

fly riiyriy tdytdij ’ 

with coefficients given by m and where p is the linear form 
of the plane containing C 002 , cqh, cioi, q is the linear form 
of the plane containing C0207 cqh, cho and r is the linear 
form of the plane containing C 200 , cioi, cuo which satisfy 

piW) = l, q{V) = l, r(C 7 ) = l, 

at the points U (intersection of the planes p,q,t), V (in¬ 
tersection of the planes p,r,t) and W (intersection of the 
planes q,r,t) given by (fS]). 

If the planes p, q, r meet at a proper point T, the quadric 
is a cone with vertex T. If T is a point at infinity, the 
quadric is a cylinder and the direction of its axis is given 
byT. 

7 . Classification of quadrics 

The classification of quadric patches is refined now that 
we know whether the quadric has a center or not: 

• A > 0 : Oval quadrics: 

— Centered: 



Figure 4: The pole of a diametral plane is the direction of the tangent 
cylinder to the quadric along its intersection. 

— Centered: One-sheeted hyperboloids. 

-- Non-centered: Hyperbolic paraboloids. 

We may tell ellipsoids from two-sheeted hyperboloids in 
other ways. For instance, if the conics at planes t,u,v,w 
are all ellipses, the quadric is an ellipsoid. We can use (HD) 
for this. 

8. Diametral planes and axes 

If a plane contains the center of the quadric, it is called 
diametral. As the center of the quadric is the pole of the 
plane at infinity, polar planes of points v at infinity are 
diametral. That is, the tangent cone to the quadric along 
its intersection with a diametral plane degenerates to a 
cylinder (see Fig.U]). The direction of the cylinder is given 
by the pole v of the diametral plane. 

We choose a basis of vectors {U,V,W}, where 

tj = Tt, V = ft, W = T#, (25) 

if U, V, W are proper points. If one of them is a point at 
infinity, we take it as vector of the basis. For instance, if 
17 is a point at infinity, we take U = U. 

The case of an improper point T can be handled simi¬ 
larly. We use in this case a basis {U,V,T} with U = Wl), 
V = hT, T = T, if W is a proper point. Otherwise, we 
choose C/ or F as origin and use W as one of the vectors 
of the basis. 

For a direction 


* Ellipsoids: detZ > 0 . 

* Two-sheeted hyperboloids: detZ < 0 . 

— Non-centered: Elliptic paraboloids. 

• A = 0 : Degenerate quadrics: 

— Cones: if the vertex T is a proper point. 

— Cylinders: if T is a point at infinity. The type of 
cylinder is determined classifying any of its conic 
sections [lfil| . 

• A < 0 : Ruled quadrics: 


V = aU -f ( 3 V 7IF, ( 26 ) 

the polar plane is a diametral plane with linear form given 

by 

_ ^ __ ^ _ 

+ (t: ?) -\- ey/3 -I- ew'l)t.{27) 

\ill/ iiy ilw J “t/ 

Before going on, we need information about the normal 
of a plane given by a linear form: 
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Lemma 1 . A plane I with linear form ap + bq + cr dt, In the case of improper T and, for instance, proper W, 
where p, q, r, t are linear forms satisfying the equations for the coordinates of iT = aU + pV + ST are 


p{W)=qiV) = r{U)=t{T) = l, 
contains the vectors 

vi = {b — evd)U + {eud — c)V, 

V2 = {a — ewd)V + {eyd — b)W, 

V3 = {ewd — a)U + (c — eud)W, 
and hence its normal vector is given by 
n= {a — ewd)U x V + (c— eud)V xW + {b — evd)W x U. 

The proof is simple, since a vector v = aU + /dV + "fW 
belongs to the plane if and only if 

0 = l{v) = (c — eud)a + {b — evd)f 3 + (a — €wd)j. 


b-eya = qi{gvuOi + gvvii + gyrS), 
c-eua = p.{guuOi + guv/d + gurd), 
d = fJ-igruO! + gTvP + grrS), 


where 

a 

b 

c 

d 


eua d- eyjd 

a fd 

) 1 


VLw 

Tlu Tty 

/ Ttw 

/3 

a 

eua d- ey/dX 

1 

Ily 


Tw J 

Tty ’ 

a 

P , 

eua d- ey/d\ 

1 

flu 

Tty 

Tw J 

Ttu ’ 


\ 5 . 


( 31 ) 


( 32 ) 


A diametral plane is called a principal plane or a plane 
of symmetry if it is orthogonal to its pole v. The principal 
axis of the quadric are the lines which are intersections of 
two principal planes. That is, the poles of the principal 
planes are the directions of the axes. 

Since for this definition we need to include a scalar prod¬ 
uct, it is necessary to provide another symmetric bilinear 
form G, such that v ■ w = G{v, w), denoting by a dot the 
scalar product. 

The matrix of such form is usually called Gram matrix 


and it is 




/ 

( guu 

guv 

guw \ 


gvu 

gw 

gyw where gAB 

\ 

V gwu 

gwv 

gww j 


in the basis {U,V,W}. 

In order to derive conditions for principal planes, we 
have to impose that the pole v be orthogonal to a basis 
of vectors of the diametral plane, which can be two of the 
ones which have been calculated in Lemma 1 , 

0 = vi ■ v = {b-eyd){guua +guv/3 + guwl) 

+ {eud-c){gyua +gyyjd + gywl), 

0 = V 2 ■ V = {a - ewd) {gyuot + gyy/d + gywl) 

+ {eyd-b){gwuoi +gwv/d + gwwl), 


These conditions can be seen as arising from an alter¬ 
native definition of principal axes as lines with direction 
given by eigenvectors of the bilinear form of the conic at 
infinity. The values of the coefficient /i are the correspond¬ 
ing eigenvalues, which are obtained by imposing that the 
system ( 1 ^^ has non-trivial solutions for a, / 3 ,7. Hence 
11 has to satisfy a cubic equation and there are in gen¬ 
eral three principal planes and axes, except for quadrics of 
revolution and spheres: 

Corollary 3 . A rational triangular quadratic patch with 
control net {coo2, coii, C020, cioi, cno, C200} and weights 
{woo2,Won,W020,wioi,wiio,W2oo}, such that the three 
boundary conics meet at a point S, which is written as in 
m, has diametral planes with linear forms given by \ 21 ^ 
with a pole v = aU + fdV + , or v = all + /dV + 7T, 

if T is a point at infinity, in a vector basis i 25 \} . defined 
with points 0) and m- 

The poles of principal planes are the directions of princi¬ 
pal axes and, if T is a proper point, they satisfy the linear 
system \ 29 i) for values of pL for which the determinant 

>~wu - - uawu Awv - nyfctv ~ r-gwv Aiv + - uaww 

^VU - ayajj - a-avu + stW “ Wvv Avw - “ uavw 

- vauu A(/v - njiiy - gauv >'UW - - vauw 


and these equations can be easily solved up to a propor¬ 
tionality factor /i, 


a — ewd 

b — Cyd 

c — eud 


p.{gwuOi -d- gwvfd + gwwl), 
fi{gvua -I- gyyfd gywl), 
fi{guuOi -I- guv/d + guwl), 


where 

a 

b 

c 

d 


1 _ 

a 


1 A 

Tlw 

Ttu 

_ 

Tty J 

Ttw 

_P__ 

a 

_ 1 3 

1 A 

Tty 


Ttw ) 


a 


_ 1 3 

1 A 

Ttu 

Tty 

Tlw J 

' Tu’ 


-X{eua + eyjd -I- ewl)- 


( 29 ) 


( 30 ) 


vanishes, with Xab ■= f-Af-sX and coefficients given by 0) 
and i 28 \) . 

If T is a point at infinity, the coordinates of the pole 
satisfy the linear system for values of p for which the 
determinant 


Auu —Tguu 
^uv — ligvu 
-ggru 


Auv — g-guv 

Ayy —ggyy 

-ggrv 


-gguT 
-ggvT 
X — ggxT 


vanishes, with 


Auu 

■ ( 


M 

, Ayy : = 



Tlu J 

Auv 

ec/ey 


ey 

1 

Q2 

“vy 


TlwTlu 



ey 

Tlw 


1 

Tly 


2 


1 

ViuXly 




















Figure 5: The vertex O of a paraboloid is the pole of a tangent plane 
with the center Z as normal. 


In general there are three different values for fi. If there 
is a double non-null solution, the quadric is a surface of 
revolution. If there is a triple solution, the surface is a 
sphere. 

The discriminant of a cubic equation a3/r^ + 02^^ + ai/i+ 
oo = 0, 


7 = 


we get the differences A = a — ewd, B = b — cyd, C = 
c — Cud for the planes which are orthogonal to the center. 
That is, we have the linear form for the tangent plane at 
the vertex of the paraboloid, except for the coefficient d. 
Since tangent planes are solutions of the implicit equation 
of the quadric in tangential form (fT 51 ) . 

0 = flu^viCeijd){Beyd) 

+ riyltwiB€vd){Aewd) 

+ Idij{A + ewd){C + eijd) 

1 

flljlty + flyflw + flw^U 
= Hi/Ity BC + Ity^lw AB + Itw^U 
+ (flijriy{euB eyC) flyllwi^V^ + ^wB) 

+ flw + ewC))d, 

the coefficients of the plane are readily obtained 


o o o o 

A = — 27 a^aQ + 1800010203 — 40301 — do^oo + O2O1, 

provides a simple way of checking whether a quadric is a 
surface of revolution: A vanishing discriminant A is equiv¬ 
alent to having a double root. 

Besides, a vanishing second derivative of the equation, 
3o3p,-|-2o2 = 0, implies a triple root and hence the quadric 
would be a sphere. 

It is easily checked that the eigenvalue p, = 0 just ap¬ 
pears for = fljjlly + Hyflw + in the case of 

proper T. That is, for paraboloids. For these surfaces the 
plane at infinity is a principal plane: The plane at infinity 
is diametral, as it comprises the center, and it is principal, 
since it fulfills (Hi). Besides, it is the polar plane of the 
center, since it is tangent to the paraboloid at the center. 

In the case of improper T, the eigenvalue fi = 0 appears 
only if A = 0 . Hence, cylinders have a null eigenvalue, 
corresponding to a pole T, the direction of the axis, which 
has no polar plane. For parabolic cylinders there is another 
null eigenvalue, since the plane at infinity is an improper 
principal plane. 

There are then just two proper principal planes for 
paraboloids, except for paraboloids of revolution. The in¬ 
tersection of these principal planes is the only proper axis 
of the paraboloid, a line with direction given by the center 
of the paraboloid. The axis meets the paraboloid at the 
center and at a single proper point named vertex. 

We may calculate the vertex solving a quadratic equa¬ 
tion, but there is a simpler way, taking into account that 
the tangent plane o at the vertex O is orthogonal to the 
center Z (see Fig. [ 51 ). We can use this property to compute 
the vertex: 

A plane with linear form ap-\-bq-\-cr-\-dt is orthogonal to 
a vector (P 51 ) if their coordinates fulfill ( 0 ^ 1 ) . Hence reading 
the coordinates of Z from (USD with, for instance, /r = 1, 

a = I}u{tly Hrv), jS = I}y{(l,w + 


VluVlyBC + Vlyidw AB + Vlwldir AC 
d = -, 

CLxj^Y (yeu B-\-evC)-\-O.YVly/ (ev-A+evv W^C/ 

a = ewd -f gwu^ui^v + ^w) + 9 wv^vifflw + ^u) 

+ gww^w{^u + ^v), 

b = eyd -b gyufluffly -b Hw) + gvv^vi^w + ^u) 

+ gvw^wi^u + ^v), 

c = end + gnij^ui^v + ^w) + guv^v i^w + ^u) 

guw^wi^u + ^v), ( 33 ) 

where a,b,c,d are given by O- 

The vertex O is just the pole of the plane o = ap-\-bq-\- 
cr -b dt: 

Corollary 4 . A rational triangular quadratic patch for a 
paraboloid with control net {coo2,coii,co20jcioi,cuo,C200} 
and weights {a;oo2 ,wqh, a;o2o,wioi,wno,a;2oo}) such that 
the three boundary conics meet at a point S, which is writ¬ 
ten as in m, has a vertex given by 

^ (bfly -b aflw)^uU -b {cfljj -b oHwjHy V 

QV — 2 e[ 7 €v'c?)“ 1 ”^V^w— 2 eY — 2 e 

{cflu -b bny)flwA^ ~ 2 d(fflij^y + ^y^w + HwHj/ )T 
H” 5 

r2[7 Q v" — 2eu €y d)-\-dY^w — 2eY d)-\-Q,w^u — 2e\y €.u d) 

where A = a — ewd, B = b — eyd, C = c — ejjd and { 23 )- 
The coefficients are given in 0 ) and il 9 \) and the points 
U,V,W,T are defined in ( 0 ) and HU) . 

Finally, the axis of a cylinder is easily determined, since 
it has the direction of T and contains the center of every 
conic section. For instance, we can use the center of the 
conic at t, which according to (HZl) is given by 

ew^W^w^emP + eyIdyLOvlXlQ2oQ + iU^U^uOJ20oR 
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C200 


Figure 6: Ellipsoid 


Figure 7: Two-sheeted hyperboloid 


9. Examples 

Now we apply our results to several quadric patches: 


Example 1. Net: 


( 0 , 0 , 0 ) ( 1 , 0 , 1 ) ( 2 , 0 , 0 ) 

(-2,2,3) (1,1,1) 


and 


weights: 


1 5 1 
1 2 
3 3 
1 


(Fig. 0); 


The normalised linear forms for the tangent planes are 
10a; + 25y — lOz 2 — x — z 12 — 3a; — 9?/ 
-32-’ = " =-16-’ 

and meet at a point (—8/7,12/7,22/7). 

The three boundary conics meet at a point, 

S = (1,0, — 1) = coo2 — Con + co2o 

coo2 — 2cioi/3 + C200 C 200 ~ 4ciio/3 + C 020 


4/3 


2/3 


and hence uju = 1, = 4/3, cOw = 2/3. 

The normalised linear form for the plane through the 
corners of the net is 


t = 


— 7y + 14z 
32 ■ 


The quadric is oval, since the coefficient A = 207/49 for 
this quadric patch is positive. 

It is not a paraboloid, since the center is the proper point 
(1/2, 2/5,—1/30). Since the boundary curves and the 
conic on t are ellipses, the quadric is an ellipsoid. One ar¬ 
rives to the same conclusion checking that detZ = 480/49 
is positive. 

The bilinear form that we get for this surface is 


Q = 


18a;^ -I- 45a;y -I- 81y^ — 81yz + 54z^ — 36a; — 90y -I- 36z 

64 ’ 

and the implicit equation, in cartesian coordinates is 

0 = 18a;^ -I- dbxy + 81y'^ — 81yz + 54z^ — 36a; — 90y + 36z. 


The three eigenvalues calculated according to Corol¬ 
lary [3] are different, y. = 1.78,0.90,0.52, and so this el¬ 
lipsoid is not a surface of revolution. 

The three principal planes are 

0.61p -f 0.379 - l-32r - 2.041, 

-0.90p -k 15.059 - 4.49r -t 5.041, 

-10.46p -k 17.249 - 0.57r - 19.121, 
with respective implicit equations in cartesian coordinates 
0.39a; -k 1.67y - 1.13z - 0.90 = 0, 

-1.32a: -k 0.72?/ -h 0.60z -f 0.39 = 0, 

-5.32a - 3.67?/ - 7.25z -f 3.88 = 0, 
and obviously meet at the center. 


Example 2. Net: 


( 0 , 0 , 0 ) ( 1 , 0 , 1 ) ( 2 , 0 , 0 ) 
(-1,1,2) (l,i,l) 
( 0 , 2 , 0 ) 


and 


weights: 


1 2 3 
1 2 
1 


(Fig.\]](: 


The normalised linear forms for the tangent planes are 
X + 3y — z a - 

P = -i-> q = — 


■ z — 2 12 — 4a — 6?/ — 5z 

- T = - 


and meet at a point (2/5, 2/5, 8/5). 

The three boundary conics meet at a point at infinity, 

S = (2,0, —4) = C 002 — 4coii-I-3 co2o 

= C002 — 2cioi + C200 = C200 ~ 4ciio + 3 co20j 

and hence = 1. 

The normalised linear form for the plane through the 
corners of the net is 

1 - 
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Figure 8: Hyperbolic paraboloid 


The quadric is oval, since the coefficient A = 3/25 for 
this quadric patch is positive. 

It is not a paraboloid, since the center is the proper point 
(7/4, —1/2, 5/2). As the boundary curves are not ellipses, 
but two hyperbolas and one parabola, the quadric is not 
an ellipsoid, but a two-sheeted hyperboloid. Accordingly, 
detZ = —16/25 is negative. 

The bilinear form that we get for this surface is 

^ + 12xy + I2y'^ + Qyz — — %x — 2dy + 8z 

48 ’ 

and the implicit equation, in cartesian coordinates is 
0 = dx^ + I2xy + 12y‘^ + 6yz — z^ — 8x — 2dy + 8z. 

The eigenvalues for the normal directions of the princi¬ 
pal planes are different, y = 0.33,-0.44,0.30, and hence 
the hyperboloid is not a surface of revolution. 

The principal planes are 

1.58p -b 0.45g - 0.91r -b 0.08t, 

0.36p -b 0.05q + O.SOr + 0.22t, 

0.79p + d.65q + 1.21r - 0.96t, 
and have the respective equations in cartesian coordinates 
0.97a; -b 1.87y -b 0.34z - 1.60 = 0, 


The three boundary conics meet at a point at infinity, 

S = ( 0 , 0 , 2 ) = Coo2 — 2 coil -b Co20 

_ C002 — 2 cioi -b C200 _ C200 ~ 2 ciio + C020 

-1/4 “ 374 ’ 

and hence uju = !> = —1/4, = 3/4. 

The normalised linear form for the plane through the 
corners of the net is 
dx — 2y — dz 


The quadric is ruled, since the coefficient A = —1 for 
this quadric patch is negative. 

It is a hyperbolic paraboloid, since the center is a point 
at infinity (0,0,—3/64), which is also the direction of the 
axis. 

The bilinear form that we get for this surface is 

^ 8z — dx^ + dy^ 

3 ’ 

and the implicit equation, in cartesian coordinates is 

0 = 2z-x‘^ + y'^. 

The eigenvalues calculated according to Corollary [3] are 
different, /i = 0,4/3, —4/3, and one of them is null, as it 
is expected for a hyperbolic paraboloid. 

The three principal planes are 

q + r + t, 

2p + q + dr + t, 

-^p + ^q + 3r + 3t. 

The first principal plane is the plane at infinity and the 
other ones have y = 0, x = 0 as implicit equations in 
cartesian coordinates. They all meet at the center. 

The vertex is the point (0,0,0), as it is clear from the 
form of the implicit equation. 


-0.05x -b 0.05?/ - 0.13z -b 0.43 = 0, 

0.75x - 0.32y - 0.39^ - 0.51 = 0, 
and meet at the center as expected. 

■ ( 0 , 0 , 0 ) ( 1 , 0 , 0 ) ( 2 , 0 , 2 ) 


Example 3. Net: 


(0,i,0) (l,i,0) 

( 0 , 1 ,-^) 


weights: 


111 
1 1 
1 


(Fig. m 


and 


The normalised linear forms for the tangent planes are 

4 — 4x -b 2z —l + 2y + 2z 
P=-z, q= -^-, r =-^-, 

and meet at a point (1,1/2, 0). 


Example 4. Net: 


weights: 


1 1 1 
1 1 
1 


(0,-l,0) (2,0,0) (0,1,0) 
(l,-ii)(l,il) 
(1,0,1) 

(Fig. 0); 


and 


The normalised linear forms for the tangent planes are 

P = l- --by, q = i- -^-y, r = x-i, 

and meet at a point at infinity (0, 0,1). 

The three boundary conics are parabolas and do not 
meet at any point. Hence the patch does not belong to a 
non-degenerate quadric. If it is a degenerate quadric, it 
is then a parabolic cylinder with direction (0,0,1). It is 
easy to check, for instance, that their respective centers 
are aligned and hence it is a degenerate quadric. 
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The normalised linear form for the plane through the 
corners of the net is 

t = —X + z. 

Using (ESI) we find a set of weights fljj = —1/2, fly = 
1 = flvy and hence the bilinear form that we get for this 
surface is 

Q = A{x + y‘^- 1), 

and the implicit equation, in cartesian coordinates is 
x = l-y‘^. 

The eigenvalues for the poles of the principal planes are 
/X = 0, 0,4, as it is expected for a parabolic cylinder. 

The principal planes are 

p + q + r, p — q- 

The first one is the plane at infinity and the second one 
is the only proper principal plane of a parabolic cylinder, 
with equation in cartesian coordinates given by j/ = 0. 


The normalised linear form for the plane through the 
corners of the net is 
^ 2 + 2x — 3z 

2 ■ 

The quadric is degenerate, since the coefficient for this 
quadric patch is A = 0. It is a cone, since the intersection 
of the tangent planes is a proper point, which is the vertex 
( 0 , 0 , 0 ).^^ 

The bilinear form that we get for this surface is 

Q = f{x^ + y^-z% 

and the implicit equation, in cartesian coordinates is 

r\ 2 I 2 2 

0 = X +y — z . 

The eigenvalues for the poles of the principal planes are 
fi = —16/3,16/3,16/3, and so the surface is a cone of 
revolution. 

The principal planes are 
p-q-r, 


Example 5. Net: 


weights: 


3 _9_ 1 

16 2 
1 


(2,0,2)(|,^,|)(i,^,l) 

(i-^4)(2>0,l) 

(Fig. \m): 


37P+ (1 - 2j)q- (1 + 7)r, 

for every value of 7 and their respective equations in carte¬ 
sian coordinates are 

z = 0, 

37 X + {2 — 'y)'/iy = 0. 


The normalised linear forms for the tangent planes are 


p = 


q = 


—X — y/iy -I- 2z —X -I- -y/Sy + 2z 


-, r = 


and meet at a point (0,0, 0). 

The three boundary conics meet at a point 


^ = (- 1 , 0 , 1 ) = 


3 C 002/8 — 9coii/8 -b C 020 


1/4 

_ 3C002/8 — 9 cioi/ 8 -b C200 

174 

and hence iZu = 1 /4 = , Wu, = 1. 


— C 200 — Clio -b C 02 O) 


10. Tensor product quadric patches 

Tensor product patches are the most common way 
to model surfaces in CAD. In particular, in some cases 
quadrics can be parametrised by biquadratic rational 
Bezier patches, 

2 2 

^ ^ {u)B] (u) 

c{u,v) = -, u,we [0,1], 

2^0 j—0 
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Figure 11: Planes on a biquadratic tensor product patch 


Figure 12: Spherical rational tensor product patch 


for a control net {cqo, cqi, C 02 , cio, cn, C 12 , C 20 , C 21 , C 22 } and 
their respective weights. 

The patch is bounded by four conic arcs with control 
polygons {coo,coi,co 2 }, { 020 , 021 , 022 }, { 000 , 010 , 020 } and 
{ 002 , 012 , 022 }, meeting two by two at the four corner ver¬ 
tices Coo, C 02 , C 20 , 022- 

Not every rational biquadratic patch is a quadric patch 
[m, but we can apply our knowledge about quadric 
triangular patches to them. 

For instance, we can take P = coo, Q = C 02 , R = C 20 
and define a triangular patch with these three corners, 
as we know that the tangent planes are defined by the 
neighbouring vertices: p contains coo,coi,cio, q contains 
co 2 ,coi,ci 2 and r contains C 20 , cio, C 21 (see Fig. E]). 

We already know the conic at u, defined by the control 
points { 000 , 001 , 002 } and their weights, and the conic at 
V, defined by { 000 , 010 , 020 } and their weights. For the 
conic at w we have the control points C 02 and C 20 and 
their weights, but we lack the intermediate control point 
Clio and the weight wno. 

In order to have a triangular quadric patch, we can use 
the other point S where the conics at u and v meet, besides 
P. We may reparametrise both conics as we did in m so 
that their weights satisfy 


Example 6. Net: 


(a,0, a) 
(2a, 0,0) 
(a,0,-a) 


with a = \l\pi and weights: 


(a, a, a) (0,a,a) 
(2a,2a,0) (0,2a,0) 

(a, a,—a) (0,a,—a) 


1 1 2 
1/2 1/2 1 
1/2 1/2 1 


fPig. [Wit : 


We use a triangular patch through P = coo, Q = C 02 , 
R = C 20 with the following control net and weights 


(a, 0,a) {a, a, a) (0,a,a) 


■ 1 1 2 ■ 

(2a, 0, 0) Clio 

5 

1/2 Who 

(a,0,-a) 


1/2 


and notice that the conic at u and the conic at v meet at 
the point S = (—•\/2/2,0, •\/2/2) and 


•S' = Coo — 2coi + 2 co 2 


Coo — Cio + C 20 12 


1/2 


There is no need to perform Mdbius transformations, 
since the weights already satisfy 

1 

Wu — oJoo ~ 2uJoi + W 02 — 1, oJv — oJoo — 2 ujiq + UJ 20 — 2) 


but the other denominator is not determined, 


^ _ ^opcpo ~ 2a;oiCoi -I- 0^02002 
OJu 


‘j-’oqCqo — 2a;ioCio -I- a;20C20 

LOy 


where denominators, if S' is a proper point, are 


aiw — aj2o — 2clIiio + UJ 02 — — 2a;iio. 

The control points and S define the planes p, q, r, t and 
their intersections. 


— oJoo ~ 2 uJoi + W02, ojy — uiQO ~ 2 uio + W20. 

Now we can dehne the plane w as the one containing 
S, C 02 , C 20 and complete Fig. [5] by computing the intersec¬ 
tion points U, V, W on plane t. 

The barycentric combinations for U, V, W provide us the 
value of uiy, and hence of wno and cuo- 

If the biquadratic patch is in fact part of a quadric sur¬ 
face, Theorem 1 provides its bilinear forms and we can 
calculate its geometric elements. We see it with an exam¬ 
ple: 


V2 V2 3 y 2 \ 



r= (1,72,72,0), 


except for W, which is a point at infinity with direction 

(1,-1,1). 

This means that Vlyy = 1, but the normalisation term 
for W 


—U! 002 aJw + 0J020‘^v + aJ 200 ^u — 2wiio — 1 
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must vanish and hence wno = 1/2, uj^ = 3/2 and the 
representative for W is the vector 


W — —LOooUJuiP + 0J02^vQ + 




The coefficients 


flu = 2, fly = 1, flw = Ij A = 2, 

yield the expression for the bilinear form for the quadric, 

^2 

Q = + — - 2pq - pr - qr + 2t^. 

The normalised forms for the planes are 


p = —^(a: + z - \/2), 
\/2 

r = —(\/2 - a; + z), t 

O 


9= ^(\/2-y-z). 


72 / 72 ^ 


and so the implicit equation for the surface in cartesian 
coordinates is 

a:^ + 7 + - 1 = 0. 


11. Conclusions 

We have derived closed formulas in terms of con¬ 
trol points and weights for several geometric elements of 
quadrics in rational Bezier form, both in triangular and 
tensor product representation. To our knowledge, these 
formulas have not been produced before. The main dif¬ 
ference with other procedures for drawing geometric infor¬ 
mation from rational triangular patches [3| is the use of 
geometric entities such as tangent planes to the quadric 
and their intersections as ingredients for obtaining bilin¬ 
ear forms, and hence, implicit equations, for the surface. 
There are many ways of implicitising a parametric surface 
, but the use of linear forms with geometrical meaning 
instead of cartesian coordinates simplifies this problem for 
quadric patches. Besides, these geometric entities appear 
naturally in the formulas for geometric elements because 
they are already present in the expressions for the bilin¬ 
ear forms for the quadric. The use of projective algebraic 
geometry allows us to perform calculations in a synthetic 
fashion, instead of resorting to cartesian coordinates. 

Additionally we classify affine quadrics using one coeffi¬ 
cient involving the weights of the patch. This can be done 
without implicitising the quadric patch [l^, [l3|; but the 
closed form for the implicit equations is what enables us 
to derive closed formulas for geometric elements. 

The results are obtained initially for Bezier triangles, 
but are also extended to quadric patches in tensor product 
form. 
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